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Abstract 

We  consider  a  single-input  nonlinear  discrete-time  system  of  the  form 
Z:  x(t+l)  *  f ( x ( t) ,u(t) } 

where  xeIRN,  ueIR,  and  f(x,u):IRN+1  -*IRN  is  a  C*  IRN- valued  function. 
Necessary  and  sufficient  conditions  for  approximate  linearizability  are 
given  for  Z.  We  also  give  a  sufficient  condition  for  local  linearizability. 
Finally,  we  present  analogous  results  for  multi-input  nonlinear  discrete¬ 
time  systems. 
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I.  Introduction  ,  . 

t  .7  ■  r  etjur  *-  .  - 

We  consider^  single-input  nonlinear  discrete-time  system  of -the  form  . 

1:  x(t+l)  =  f ( x ( t ) >u( t) )  (1) 

where  xeIRN,  ueIR,  and  f(x,u):  IRN+1  ■+1RN  is  a  C*  IRN-valued  function. 

Many  authors  have  studied  (local  or  global)  linearization  (Cheng  et.  al. 

1985,  Hunt  and  Su  1981,  Jakubczyk  and  Respondek  1980,  Krener  1973,  Su  1982) 

and  approximate  linearization  (Krener  1984)  by  state  feedback  and  coordinate 

change  for  nonlinear  continuous-time  systems.  In  this  paper  .we  discussei. 

necessary  conditions  and  sufficient  conditions  for  local  linearization  and 

approximate  linearization  by  state  feedback  and  coordinate  change  for  nonlinear 

discrete-time  systems.  Other  related  work  on  nonlinear  discrete-time  systems 

can  be  found  in  (Grizzle  1985a,  1985b,  Grizzle  and  Nijmeijer  1985,  Monaco  and 

Normand-Cyrot  1983a,  1983b).^W^'*V-  .  gj,\  ,.>/!  )  . 

Definition  1:  A  point  (x_,ua)  such  that  f(xa,u.)  =  xa  is  called  an  equilibrium 

“  e  e  e  e 

point. 

Now  consider  the  following  linear  discrete-time  system  Zq. 
lQ:  y(t+l)  *  Ay(t)+bv(t)  =  g(y(t) ,v( t) ) , 


where 
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(NxN  matrix) 


(N*l  matrix) 


Similar  to  the  continuous-time  case  (Krener  1984,  Su  1982),  we  can 
define  local  linearlzablllty  and  approximate  llnearizabllity  for  a  discrete- 
time  system.  Let  (xg,ue)  be  an  equilibrium  point  of  Z. 


i 


□ 

□ 


'odes 


Definition  2:  I  is  said  to  be  locally  linearizable  at  (xe«ue)  if  there 
exist  an  open  neighborhood  U(<=IRN+1)  of  (*e,ue)  and  a  diffeomorphism 
T:U-*-T(U)  such  that  (i)  T  =  (T^  ,T2 . T^)  are  functions  of  Xj,x2,...,xN 

only,  (ii)  T(xe»ue^  S°(N+1  )x1*  and  ^ii1^  ^f  =  g«>T. 

If  we  let  (y(t)Tv(t))T  =  T(x(t) ,u( t) ) ,  then  y(t)  and  v(t)  satisfy  ZQ. 
Definition  2  indicates  that  we  want  to  find  a  diffeomorphism  T  such  that 
the  diagram  in  Figure  1  commutes.  Once  we  find  such  a  diffeomorphism, 
we  can  apply  linear  system  theory  instead  of  nonlinear  system  theory. 
Definition  3:  Z  is  said  to  be  approximately  linearizable  with  order  p  if 
there  exist  an  open  neighborhood  U(c=  DT  )  of  (xg,ue)  and  a  diffeomorphism 
T : U  -*T(U)  such  that  (i)  T  *  (T^  ,T2,. . .  ,TN)  are  functions  of  x-j  ,x2,. . .  ,xN 
only,  (ii)  T(xe*ue)  *0(N+])xi  •  and  (HI)  T-f  =  g»T  +0(x-xe,u-ue)p+1 . 

Thus,  in  Definition  3  we  consider  the  following  nearly  linear  discrete¬ 
time  system: 

Zq:  y(t+l)  =  Ay ( t )  +bv(t)  +0(x-xe>u-ue)p+^ 

where  the  N*N  matrix  A  and  N*1  matrix  b  are  the  same  as  Zq.  Clearly,  local 
linearizability  at  (xe,ug)  implies  approximate  linearizability  with  arbitrary 
order. 

In  Section  II,  some  background  material  is  reviewed  and  notations  are 


defined.  In  Section  III,  necessary  and  sufficient  conditions  for  approximate 
linearizability  will  be  given  for  the  system  (1).  Also,  we  will  give  a 
sufficient  condition  for  local  linearizability.  We  can  define  local 
linearizability  and  approximate  linearizability  for  multi-input  discrete-time 
systems  similarly  to  Definitions  2  and  3.  Then  the  multi-input  case  will  be 
discussed  in  Section  IV. 


II.  Preliminaries 

In  this  section,  notations  and  definitions  to  be  used  later  will  be 
mentioned.  The  Kronecker  product  Is  very  useful  in  the  field  of  matrix 
calculus  (Graham  1981).  First,  define  the  Kronecker  product®  by 


A  ®  8  = 
p*q  m*n 
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where  a^  is  the  (i  ,j)-component  of  the  pxq  matrix  A. 

Define  the  derivative  of  a  matrix  with  respect  to  a  matrix  by 
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We  also  define 


Djj B-  B,  D^B  =  DaB,  and 
Da+1  B  =  Da(D^B)  for  1  >1. 

N  k 

Let  h(x)  be  a  scalar  real  valued  function  of  xeIR  .  Then  (Dx  h)(x)  and 

(D  *L  h) (x)  are  Nkxi  and  lxNk  vectors,  respectively. 
x1 

Fact  (Vetter  1970,1971):  Using  the  definition  of  Kronecker  product  and 
derivative  operations  on  matrices,  Taylor's  formula  can  be  expressed  by 

t  ,  k 

h(x)  *  h (0)  +  l  j^y  (D  kt  h(x))xsQ(x®x®...  ®x) +R£+1(x*), 
k=l  ’  x 


where  R«a1(x*)  is  a  remainder  term. 


k  k 

Now  define  the  N  *N  permutation  matrix  U.  s  .  as  follows: 

1 1 » ■  2* 

the  (a.  -l)Nk-1  +  (a.  -1 )Nk~2  + . . .  + (a.  -1 )N  +  a .  -th  column  of  U4  . 

M  ’2  1k-l  ’k  V’2 . \ 

is  the  (a1-l)Nk”^ +(a2-l)Nk"2 +  ... +(ak_^-l)N +ak-th  column  of  NkxNk 

identity  matrix  (I  l  u),  for  1  < a, ,a0,. . . ,a.  < N  (the  (a.)  are  related  to  the 
N  xN  -  1  2  k-  1 

"base  N"  representation  of  the  column).  Here  (i i  ,i2,. • • *i k )  is  a  permutation 
of  {l,2,...»k}.  For  example,  when  Ns2  and  k=3, 

U1 23  =  *8x8 
and  r 

1  0  0  0  0  0  0  0 

0  0  0  0  1  0  0  0 

0  0  1  0  0  0  0  0 

0  0  0  0  0  0  1  0 

U321  ”  01000000 

0  0  0  0  0  1  0  0 

0  0  0  1  0  0  0  0 

0  0  0  0  0  0  0  1 


Let  A  be  a  pxN  matrix.  Define  the  operator  by 

®  A  =  A  (  ^i  i  ^  ) 

k-  all  permutations  1l*12”*'*1k 

{i  1 ,1 2 » •  •  •  »1  k ^  {1*2*. ..»k} 

3 

For  example,  when  A  is  a  pxN  matrix, 

3*j  A  *  A^U123+U132+U213+U231  +U312+U321^ 

Let 

//3[f  \  /  3f  \  /  3f  \  (ilf.\N“l  f  — 1 

V  3u  '(0,0)’  l  3x  '(0,0)  1  3u  '(0,0) . 1  3x  '(0,0)  '  3u  '(0,0) 

|j 

be  linearly  independent;  that  is,  they  form  a  basis  for  IR  .  Define 
C:IRN-*’  IR  by  c(v)  «aN,  where  v  Is  a  lxN  row  vector  and 
WT  ,  r  „  /  3f  \1-1  /  3f  i 

r  3x)  o,o)1  auMo.o)1 


III.  Single-Input  Case 

In  this  section,  our  main  results  will  be  given.  If  f(x,u)  has  an 
equilibrium  point,  without  loss  of  generality,  we  can  assume  that  f (0 ,0 )  =0; 
for,  if  not,  let  x  =  x-xg  and  u  =  u-ufi.  Then  x(t+l)  * f(x(t),u(t)) 

*  f(x+xe,u+ue)-xe  with  f(0,0)  *0. 

Let 

^(x.u)  * f (x,u) , 

f1+1(x,u)  =  f(f1(x,u),0),  for  1<1<N-1. 

f^x.u)  represents  the  effect  of  an  input  u  at  t*0  on  the  state  at  t=i.  ffx.u) 
is  essential  for  solving  many  problems  arising  in  discrete  time  nonlinear  systems. 
Lemma  1 :  Z  is  locally  linearizable  at  (0,0)  if  and  only  if  there  exists  a 

oo  N 

C  function  h:W(<=  DT)-*  IR  such  that  (i)  W  is  an  open  neighborhood  of 
OeIRN,  (ii)  Du(h<>f^)  =0  on  some  neighborhood  of  OeIRN+1 ,  for  1  £i  _<N-1, 


(111) 


3x'x=0 


(1^)(0,0)  *  °* 


«•  1  3x  ;(0,0H 

(iv)  (Du(h«f^))(o,o)  t  0,  and  (v)  h(0)  *0. 


Proof:  Necessity:  Suppose  that  Z  is  locally  linearizable.  Thus  we  have 
a  diffeomorphism  T.  Let  h(x)  «T^(x).  (Since  T^(x,u)  depends  only  on  x, 
we  can  write  T^(x)  instead  of  T^(x,u).)  Note  that  Tg'T^f.  Since 
Du^T2^  on  some  "^hborhood  of  the  origin,  Du(Ti ®f )  =0  on  some 
neighborhood  of  the  origin.  From  now  on,  for  convenience,  we  will  omit 

"on  some  neighborhood  of  the  origin".  Note  that  T^  “Tg^f  *T^»f  . 

2  a2 

(Actually,  we  can  write  Tj*Tj*f  ,  because  T^»f  depends  only  on  x.  But  f 

Ao 

is  used,  for  consistency  of  notation.)  Since  D^T^)  =0,  Du(T^»f  )  =0. 
Proceeding  in  this  manner,  since  T^  *TN _^«f  »  . . .  «T^#fN"^  and  Dy(TN)  =0, 
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Du(Tl.f  "  )  =0.  Thus  we  have  shown  that  Du(T1ofl)  =0,  for  1  <_i  <_N-1 . 
Since  T  is  a  diffeomorphism,  T.  -T^?1’1  for  2  <  i  <  N+l ,  and  T]  .Tg.. . .  ,T^ 
depend  only  on  x, 


(  — ) 

1  3x'x=0 


'  3x  '(0,0) 


1  3x  '(0,0) 


and  Du(h®fN)  to.  Since  ^(0,0)  =0,  h(0)  =0. 

Sufficiency:  Suppose  that  there  exists  h : IR  -*-IR  satisfying  the  given 
conditions.  Let  T.(x)  =  hof1-^,  for  1  £i  £N+1.  Then  it  can  be  easily 
checked  that  T«f*g°T  and  T(0,0)  *0.  Since  det((  3^lu)  0)^0, 

there  exists  an  open  neighborhood  U  of  (0,0)  such  that  T:U-*-T(U)  is  a 


diffeomorphism  by  the  inverse  function  theorem. 


(Q.E.D. ) 


Let  5  •  (  *)• 
Lemma  2:  I  is  appi 


linearizable  with  order  p  if  and  only  if  there 


exists  a  C°°  function  h:W(c  BT)  -*IR  such  that  (i)  W  is  an  open  neighborhood 
of  0eIRN,  (ii)  (D^D^h®?1 ) ) (0 ,o)  “ °(N+1  )0xl  for  1  -1  -N_1 


(iii) 


(  — ) 
v  3x 'x=0 

'  Sx  '(0,0) 


(  ath-f1*'1)  , 

'  3x  '(0,0)* 


(iv)  (0u(h^N)){0>0)  1*0,  and  (v)  h(0) -0. 


Proof:  Necessity:  Suppose  I  is  approximately  linearizable  with  order  p. 


Let  h(x)  =Tj(x).  By  definition,  T^ffx.u)  =T2(x) +0(x,u)p+1.  So 

(D^(Du(h°f))(0j0)  =0  forO<j<p-l.  Note  that  T^?2  =T2»f  +  0(x,u)p+1 . 

Since  T2»f(x,u)  *=T3(x) +0(x,u)p+1  by  definition,  T^f2  =T3(x) +0(x,u)p+1. 

Thus  (D^Du(h°f  ))(g,o)  f°r  Proceeding  in  this  manner,  we 

can  show  that <D^Du( h**? 1  ))(0>0)  =0  for  1  <i  <N-1  and  0  <j  <p-l.  Note  that 
3T. 

(  lx"  ^(0,0)  *  ^  3x  n\  forl<i<N.  Since 


(0,0) 


(  —L  ) 

1  3x  '(0,0) 
3T2 

*  lx  ho,o ) 


(  &  ) 

v  3x  ;x=0 

/  3(h°f)  » 

K  3x  '(0,0) 


{  !!n  i 
v  3x  '(0,0) 


( ath^'1)) 

1  3x  ' 


(0,0) 


It  can  be  easily  shown  that  T^f^x.u)  =  TN+1(x,u)  +0(x,u)p+1.  Thus 

(Du<^N))(0,0)  *(dutn+i(x’u))(o,o)  Ffna11**  »»(o)  =  V0)  =0- 
Sufficiency:  (by  construction)  Let 


T-|(x)  =  h(x) 


1 


k  times 


T2(x)  ■  J  (DxT(h®f ) )(Q,0)  (x®x®***®x) 


T3(x)  '  J,  FT  (D^T(h’^))(0,0)  (x®x®---®*) 


tn(x>  ■  J,  FT  (DxT(h,fN’1))(0,0)  <****•••  ®*> 


Tn+,(x,u)  =  T,.f"(x.u). 


Then  it  can  be  easily  checked  that  T  as  defined  above  satisfies  the 
conditions  of  Definition  3. 


Now  note  that 


(D^Du(h°f  )){0i0j  =  (Bm,£^(0,0)(Dx  1  h)x=0’ 


where 


C  *  dj(du?1it 


,•  J1  2iT  kl‘k2,  -iT  k?"k3 

C"  k[.|  k^,"k£l,  °5  (V  «°t  <¥  »Df 

•(Dj-f'7*...  eD^'1  k^(D5flTgD^  Duf,T)  ...)),  for  l<£<m. 
(For  a  proof  of  (2),  see  the  Appendix.)  Let 


A11  A12  Alk 

A21  A22  ...  A2k 


Akl  Ak2  Akk 


where  A.  ,  is  an  ( N-l  )(N+1  j’xN^  submatrix  defined  by 

*  J 

Aij  =0(N-l)(N+l)ixNJ+l  ,  lf  1  <J» 

{BiV(0,0) 

A..  =  •  (Bij^(0,0)  #  ifi>j. 

J  ( j+1 ) ! 

Let  the  (N-l )(N+1  *1  vector  31  be  defined  by 


■w 


(B1 •0)(0.0> 


L(Bi,0^(0,0) 


) 


(See  Section  II  for  the  definitions  of  ®  and  ;.)  Also,  let  Bk  = 

1T  2T  kTT  j‘ 

(6  6Z  ...  6  )  . 

With  these  preliminaries,  we  can  state  our  main  theorems. 

Theorem  3:  Z  is  approximately  linearizable  with  order  p(>.  2)  if  and  only  if 

/  •  x  I  /  3f  »  /  3f  v  /  3f  ,  /  3f  /  3f  »  | 

m  V  3u  '(0,0)’  [  3x  '(0,0)'  3u  '(0,0) . '  3x  ' (0,0)'  3u  '(0,0) /  are 

linearly  independent,  and 

(ii )  Bp_1  €  Image (Ap  l). 

Proof: 

Necessity:  Suppose  that  Z  is  approximately  linearizable  with  order  p. 

Then  there  exists  a  function  h(x)  satisfying  (i)-(v)  of  Lemma  2;  in 

;i 


particular,  (Du(h°f  ) )(0>0)  s  <  3x  *x=0  ^  3x  ^(o!o)  *  3u  ^(0,0)  =  0  for 
1  <i  <N-1,  and 

Assume  that  {(|J  )(0,o)^lx  ^(0,0)^lu  ^(0,0) . {fx  ^(O.O^lu  *(0,0)} 


are  not  linearly  independent.  Then  there  exists  k  such  that  1  £k£N-l  and 

,k-l 


^  3x  ^(0,0)  ^ 3u^(0,0)  "  aj^3x^(0,0)^3u^(0,0)  for  some  constants  {otj}j=0‘ 

Thus  <f£)N  ^(q.o)  *  “j^lx^olo)  ^fu*(0,0)’  and 


n  (— /  3f  \ 


^3x  x=0^9x^(0,0)^3u^(0,0)  =  °*  This  is  a  contradiction,  which  implies  that 

'■'3u'(0,0)  *'3x'(0,0) 'Sir  (0,0)  * "  ‘  *  3x'(0,0)  3ii)(0,0)}  are  linear^ 
Independent. 

Recall  that  (Dg(h°^  >>(0,0)  =  (^(h*?1  ))(0>0)  s  0  for  m=0  and 

l<i<N-l.  Since  {^)(0,0)’(l7)(0,0)(lu)(0,0)*-*-’(l7)(o!o)(l^)(0,0)1 
are  linearly  independent,  (Dxh)x=0  is  uniquely  determined  up  to  a  constant 
multiple  (i.e.,  ( Dxh ) x_q  =  ac,  where  the  scalar  a(^0)  is  arbitrary  and  the 

Nx)  column  vector  c  satisfies  <=T(|£) >Q) (|J) (Q >Q)  =  0.  for  0<i  <N-2, 

and  c^(—  1^”^  1  -  11 

and  c  '3x'(0,0)'3u'(0,0)  ~  1}* 

Now  by  (ii)  and  (iv)  of  Lemma  2,  (D^DJh.f1  ))^Q  0)  88  0  for  1  £m  £  p-1 
and  1  <i  £N-1.  From  (2),  we  obtain 


<Dxh>x*0 


<Dxh>x=0 


Rh)> 


<°K-o 


P-1,1 


P-1 ,2 


p-1 ,p-l 


>1,0 


RN-1  ,  y 
Bp-1 ,0  (0,0) 


Since  ^B1 ,0^(0,0)^Dxh^x=0  =  ^DCDu^T^(0,0)^Dxh)x*0  =  a  ((DCDu?T)(0,0)^ 


<°H-0 


(0,0) 


Thus  the  right-hand  side  of  equation  (3)  is  -cxC (61  )T,. . .  »(BP_1  )T]T  “-ae^ . 
It  follows  that  6  ,  Is  in  the  Image  of  the  matrix  on  the  left-hand  side 

p-i 

L 

of  (3).  However,  the  {Dx  h}  are  constrained  because,  for  example, 

2  2  2 

3  h/3x. 3x.  *  3  h/3x.3x. .  Hence,  the  stronger  condition  @  ,  e Image  (A  ,) 

i  j  j  1  p— i  p—  i 

holds,  as  is  proved  in  the  Appendix  in  Lemma  A. 2. 

Sufficiency:  Suppose  that  (i)  and  (11)  above  are  true.  By  (1),  there 

exists  an  Nxl  vector  C1  such  that  c{(|^)(o,o)^lu^(0  0)  *°  for 

and  Cl(lxJ(o!o)^(0,0)  By  there  ex1st  C2*C3 . Cp  such  that 


2!  C2 
IT  C3 


pTCp 


where  is  an  N^xl  vector.  Let 


h(x)  *  J  tt  cT  (x®x®  ...  ® x) 

i«i  1!  1  —  -  —  ■  — 

1  1  1  times 


Then  It  can  be  easily  checked  that  (D^Du(h»f^ ))(Qto)  *®»  *or  ^ 

,xd  0<J  <o-l.  Clearly  <Du(h.fN))(M)  •  <oK-0<!j><o!o)(!j><0.0)  f0' 


Now  assume  that 


det 


(3h) 

v3x'xsO 


<^>(0.0) 

<^(0,0, 


=  0. 


k-1 


Then  there  exists  k  such  that  1  _< k  <_N-1  and  (°xh^x«0 ^1x^(0  0)  =  ^ 
■“lKh,x=0(!x)(0,0)’  for  some  (Vl-0'  Thus  (Dxh*x»0(lx,(olo)(lu)(0,0) 

*  V  «,(D>)I.o'f>(o!o)+,(f)(o,o)  *0'  This  1s  a  “ntr»<:tion-  "tl,ch 

implies  (iii)  of  Lemma  2.  Hence,  by  Lemma  2,  Z  is  approximately 
linearizable  with  order  p.  (Q.E.D.) 

Remark:  Z  is  approximately  linearizable  with  order  1  If  and  only  If  (1)  of 
Theorem  3  holds,  just  as  In  the  continuous  case  (Krener  1984). 

Now  a  sufficient  condition  for  local  linearizability  Is  given  In  the 
following  theorem. 

Theorem  4:  Suppose  that  f(x,u)  of  Z  is  an  analytic  IR  -valued  function. 

Z  is  locally  linearizable  at  (0,0)  if 

W  ^fu^O.Oj’^lx^O.O^lu^O.O) . <U)(01,0)(l^>(0,0)}  are  linearly 

independent,  and 

( 1  i )  there  exists  k(<«)  such  that  8^espan(C^)  for  all  £>1,  where  c£  Is 
composed  of  the  first  k  columns  of  A^. 

Proof:  By  (1),  there  exists  an  N*1  vector  c-j  such  that  c^  (^(o.o^Iu^O.O) 


T ,3f \N-1  ,3f \ 


0  for  0  <_1  <N-2  and  c,' (f^;o)(£)(0j0) 
C2*Cj,...,Cj  such  that  j  <»  and 


1.  By  (11),  there  exist 


14 


»a'l 


I  Ml 
I '  < 


I®1 


U 1 


& 


nil 


8 


IS 


vl« 

\v 


i 


m 


a 


m 


?i 


_LC 

2!  c2 


J_c 

31  c3 


•Bt»  for  l>  j. 


j!  Cj 


where  c.  Is  an  N  *1  vector.  Let 


•  i  J-cT 

ill  1 


1  terms 


(X9XS...  «x). 


Then  it  can  be  easily  checked  that  ( D^Du(h°?i ) ) (o ,0)  =0  ^or  1  — *  and 
s>0.  Since  both  h(x)  and  f1  are  analytic,  h-f^x.u)  are  analytic,  for 
1  <n-1.  Thus  D^h*?1)  =0,  for  1  <1  <N-1.  As  in  the  sufficiency  proof 

in  Theorem  3,  it  can  be  shown  that  h  satisfies  the  other  conditions  of 


Lemma  1 . 


(Q.E.D. ) 


It  is  easy  to  see  that  (ii)  of  Theorem  4  implies  (ii)  of  Theorem  3. 


Example  1;  Consider  the  following  discrete-time  nonlinear  system 


x-j  ( t+1 ) 


x2(t+l) 


x2(t)+2x1(t)u(t)+u(t)< 

XjCO+uft) 


f ( x ( t ) ,u(t) ) 


s,"“  (f)(0.0)  ■  <  1 )  *nd  <£>(0,0)<&><0,0)  *  <  J)-  <f)  °f  Theorem  4  ,s 


satisfied.  Note  that 


»’■  «wT>(o.o))  ■  '([1  °J)  ■  [\] 


and  B1  *  0  4  for  i  >2. 


4fe 

$ 

AS 


a 


m 


% 


*Bn*(o,o) 


'  (VT)(0,0)  ®  <D/>(0.0) 


0  T 

10  ®  (0  1) 

0  1. 


0  0  0  1' 
0  10  0 
0  0  0  1 


t  ro  o  o  2"| 

‘i-  2*  (B'i)o-°) ■  Lsi^iJ 


It  can  be  easily  checked  that  all  elements  of  the  4-th  column  of 


are  0  for  i  >2.  Thus  B^espanCC^)  for  £>.1,  where  C^  is  composed  of  the 


first  4  columns  of  A^.  Therefore  (11)  of  Theorem  4  is  also  satisfied. 
Hence  Z  is  locally  llnearlzable  at  (0,0).  Actually  we  can  construct  a 
diffeomorphism  T  = (T^ .Tg.T^)  in  the  way  that  is  given  in  the  proof  of 


Theorem  4.  Since 


0 

_L  o 
2!  0 
-2 


_L0 

3!  U8*l 


*  -8^  for  l  , 


L  TFTJT  V+ixl  J 

cJ2  =  (0  0  0  -2).  Clearly  c{  -  (1  0).  Thus 


2  2  2 
Tl(x)  =  X1  ■  2l  x2  *  xrx2 


T2(x)  =  T^ffx.u)  *  x?-x2 


T3(x,u)*T.|»f2(x,u)  *  x1  +u-(x2+2x1u+u2)2. 


Example  2:  Consider 


x^t+l)]  f  x2(t)+2x1(t)u(t)+x1(t)2u(t)+u(t)< 


x2(t+l) 


x] ( t)+u( t) 


f ( x ( t ) »u(t) ) 


Clearly  (i)  of  Theorem  3  Is  satisfied,  because  (|J)(o  o)  and 
(lx,(0,0)(lu)(0.0)  are  the  same  as  in  ExamPle  Since  <D5DufT)( 0,0) 


[o  o].  61  •  C((D£o/)(0j0))  -  [o]  .  Since  (0 


!DufT)(0,0) 


0  0,  B2  =  (2  0000000  0)T. 


o  0  0  2  1 

\n  •  0  1  1  0  I 

1  Lo  0  0  2J 


tB22*(0,0)  *  (°CfT)(0,0)  ®  (0CfT)(0,0)  ®  (DufT)(0,0) 


[!  ?]•[!  W 


(0  1) 


0  0  0  0  0  0  0  1 

0  0  0  0  0  1  0  0 

0  0  0  0  0  0  0  1 

0  0  0  1  0  0  0  0 

0  1  0  0  0  0  0  0 

0  0  0  1  0  0  0  0 

0  0  0  0  0  0  0  1 

0  0  0  0  0  1  0  0 

0  0  0  0  0  0  0  1 

• 

0  0  0  0  0  0  0  6 

00020220 
00000006 
00020220 
02202000 
00020220 
00000006 
00020220 

.0  0000006 


•’  A22  '  3,  (B22J(0*0) 


(B2!*(0.0) 


<  WT  *  °/»(o.o) 


+  (VT®WT)(o,o) 


<D5fT®°Ao,0)  + 


V  *  Dx10ufT'1 
vT  *  VufT 

LVT  *  W*  . 


+  (D?f' 


DCDufT)(0,0) 


0 

2 

0 

2 

0 

2 

0 

2 

0 


0 

0 

2 

0 

0 

0 

2 

0 

2 


4 

0 

4 

0 

0 

0 

4 

0 

4 


'21 


*  (B21 }(0,0) 


0  1 

0 
0 
0 
0 
0 
0 
0 
OJ 


0 

4 

0 

4 

0 

4 

0 

4 

0 


4 

0 

6 

0 

0 

0 

6 

0 

6 


4 

0 

6 

0 

0 

0 

6 

0 

6 


0 

0 

0 

0 

0 

0 

0 

0 

0 


Since  B-j  €  Image  (A^),  I  is  approximately  linearizable  with  p=2.  However, 
since  B2  ^  Image  (A2),  It  is  not  approximately  linearizable  with  p=3.  Thus 
it  is  also  not  locally  linearizable.  Let 

2 


1 


xl‘x2 


x2  ~X1 


T3  »  Xj  +u  -  (x2  +2x^u  +x.|2u  +  u2)2 


Then 


<11 


y,(t+l)  « 


T1(x(t+l))  =  x2(t)  +2x1(t)u(t)  +Xl(t)2u(t)  +  u(t)2  -  (x,(t)  +u(t))2 
*  x2(t)  -x^t)2  +x1(t)2u(t) 


=  T2(x(t))  +0(x,u}3 
a  y2(t)  +o(x,u)3 

T2(x(t+1))  ■  T3(x( t) ,u( t) )  ■  v( t) 


IV.  Multi-Input  Case 

The  results  in  Section  III  can  be  easily  generalized  to  the  multi¬ 
input  case.  Thus,  in  this  section  we  give  (without  proof)  a  sufficient 
condition  for  local  linearizability  and  a  necessary  and  sufficient 
condition  for  approximate  linearizability  by  state  feedback  and  coordinate 
change  for  a  multi-input  nonlinear  discrete-time  system  (for  proof,  see 
Lee  1986). 

Consider  a  multi-input  nonlinear  discrete-time  system  of  the  form 

Z:  x( t+1 )  =  f (x(t) »u( t) )  (< 

where  x(t)  eIRN,  u(t)  elR"1,  and  f(x,u):  IRN+m  -*>IRN  is  a  C°°  IRN- valued 
function.  Also,  consider  the  following  multi-input  linear  discrete-time 
system 

2q:  y(t+l)  *  Ay(t)  +Bv(t)  *  g(y(t),v(t)), 

where  y(t)  elRN,  v(t)  ellf,  A  *  block  dlag  {A^  .A22,. . .  .A^} , 


0  10.. 
0  0  1.. 


0  0  0  .. 

0  0  0  .. 


(K.xKj  matrix). 


•  M 

,1, 

B  «  block  dlag  {b.j,b2,...,bm},  and 

bj  *  (0  ...  0  1)^.  (K^xl  matrix) 

Definition  4;  l  Is  said  to  be  locally  llnearizable  at  (xe,ufi)  If  there 
exist  indices  an  open  neighborhood  U(<=  IRN+m)  of  an  equilibrium 


S 

\wr  ' 

.1  V. 


point  (xg,ue)  and  a  diffeomorphism  T:  U-*-T(U)  such  that 


(i)  T  =  (T.|  ,T2.. . .  ,Tn)  are  functions  of  x.|  ,x2,. . .  .x^  only, 

(ii)  T(xa,ua)  =  0lMimWl,  and 


(N+m)*l 


(iii)  T°f  =  g»T. 


If  we  let  =  T(x(t) ,u(t)),  then  y ( t)  and  v(t)  satisfy  the 


relation  E, 


Definition  5:  E  is  said  to  be  approximately  linearizable  with  order  p 


if  there  exist  indices  an  open  neighborhood  U  (c=lRN+m)  of  an 


equilibrium  point  (x  »u  )  and  a  diffeomorphism  T:  U-*-T(U)  such  that 


(i)  T  =  (T-j  ,T2.. . .  ,Tn)  are  functions  of  x-|  ,x2,. . .  ,xN  only, 

(ii)  T(x..uJ  -  0(N+m)xi,  and 


,  and 


(iii)  T»f  *  g«T  +0(x-xe,u-ufi)p  . 


Thus  in  Definition  5  we  consider  the  following  nearly  linear  multi - 


input  discrete- time  system: 


y(t+l)  =  Ay(t) +Bv(t) +0(x-xe,u-ue)f 


where  the  N*N  matrix  A  and  N*m  matrix  B  are  as  in  EQ. 


Now  we  state  the  generalized  version  of  Lemmas  1  and  2  and  Theorems  3 


and  4.  Just  as  in  the  single-input  case,  we  can  assume  f(0,0)  =0  without 


loss  of  generality,  if  f  has  an  equilibrium  point.  Also,  define  fix.u) 


in  the  same  way  as  in  the  previous  section. 


Lemma  5:  E  is  locally  linearizable  at  (0,0)  if  and  only  if  there  exist 
tKi}i*«i  and  c*  functions  hj(x),h2(x),...,hm(x):  W(<=IRN)-»  IR  such  that 


(i)  W  is  an  open  neighborhood  of  OeDT, 


(ii)  D,,(h . »f 1 )  =0,  for  1  <j  <m  and  1  <1  <k,-l. 


and 

(v)  h.(0)  =0  for  1  <  j  <m. 

J 

Let  £  =  (xT,uT)T.  Thus  £  is  a  (N+m)*l  vector. 

Lemma  6:  Z  is  approximately  linearizable  with  order  p  if  and  only  if 
there  exist  and  C*  functions  h1(x),h2(x),...,hm(x):  W(c=IRN)-*-  IR 

such  that 

(1)  W  is  an  open  neighborhood  of  OeIRN, 

(ii)  (DjrDu(hj*^))(o  o)  *or  1  — J  ^  and  and 

(iii) ,  (iv),  and  (v)  of  Leirma  5  are  satisfied. 

Let 
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E=  ^311^(0,0)*  (3x)(0,0}(3u1)(0,0)’---’(3x)(0,0)(3u1)(0,0)’  ^u^O.O)’ 

,  3f ,  /_3f_\  /SfxSn'1  ,  3f,  , 

•'•^3x;(0,0)^3u2;(0,0)**"a3um;(0,0)’"-^3xJ(0,0)'3um;(0,0)i 

Ei  ={(3u^(0,0)*  ^(O.OJ^CO.O) . (U)(0,0)(^J7)(0.0)’  ^(O.O)’ 

•••’(lr)(0,0)(^)(0,0)’---’{^;)(0,0) . (l7)(0,0){^f)(0,0)}’ 

l  ™  1 ,  • » .  ,m . 


Suppose  that  the  elements  of  E  are  linearly  independent;  that  is, 

N  1 

they  form  a  basis  for  IR  .  Let  a.  =  l  K,-  for  1  <i  <m.  Define 

j=l  J 

C1  ( v ) :  IRN  -*-IR  by  s’(v)  =  a  ,  where  v  is  a  1*N  row  vector  and 

°i 

T  _  ,  3f ,  .  .  /3f\Kr1  /3fv  .  /  3f , 

v  Vsu^fo.o)  **•  “o^ax^o.oj^u^o.o)  “cyi^u^o.o)  ••• 

+  °o2(lx)(0i0)^)(010)  +*”  +\.1+1(^)(0*0)  +  ”’ 

.  /  3f  \  *Sn  /  3f  \ 

+  aom'3x'(0,0)'3um  (0,0)' 

Also,  define  ^ ;  IR^+IRP,  for  i=l ,2,.. . ,m,  by 


^(V)  =  (C1  (v, )  C1(v2)  ...  cVp))T 


where  v.  is  the  j  row  of  V.  Let 

J 

■(°KfT,(o.o) 

„  i  .  (D^D.f2T)/n  \ 


FT  VI  v  I  T?  V  t^rrwr- 


»  vwirw  V*  v*»  ^  w  ^  »  "’I  ’ 


s 

;V 

r-' 

,% 

ii. 


Also,  let  = 

((sj)T 

(sj)T  ... 

(<)T 

r  d1 

un 

u12 

D1  1 
ulk 

U21 

1*22  •  •  • 

• 

D1 

u2k 

Di 

L  kl 

uk2  •  •  * 

Dkk  . 

Let 


where  .  =  0  .  ...  if  i  <j,  and 

1J  m(k£-1)(N+m)1xNJ  1 


(Bi  j^O.O) 


of .  =  • 

1J  (j+D! 


(Bi  j^O.O) 


Bi j  J(0,0)j 


if  i  >j 


Theorem  7:  £  is  approximately  linearizable  with  order  p(  >_2)  if  and  only 
if  there  exist  (K . such  that 

(i)  the  elements  of  E  are  linearly  independent, 

(ii)  span  E..  *  span  (E^fl  E)  for  1  <1  <.m,  and 

(iii)  Yp_-|  e  Image  ( 0^_  1 )  for  1  <£<m. 

Remark:  £  is  approximately  linearizable  with  order  1  if  and  only  if  (i) 
and  (ii)  of  Theorem  7  hold,  just  as  in  the  continuous  time  case  [6].  If 
m=l  (single-input  case),  K^*N.  Thus  (i)  of  Theorem  7  is  the  same  as  (i) 
of  Theorem  3.  Since  E^  =E^fl  E,  (ii)  of  Theorem  7  is  trivially  satisfied. 
Since  the  operator  ^  is  the  same  as  the  operator  c  in  the  previous  section, 
y1  .  =  8  Since  D1  ,  *  A  . ,  (iii)  of  Theorem  7  is  the  same  as  (ii)  of 

'p-1  p-1  p-1  p-1 

Theorem  3.  Therefore,  Theorem  7  is  a  generalized  version  of  Theorem  3. 

Now  a  sufficient  condition  for  local  linearizability  is  given  in  the 
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following  theorem. 

Theorem  8:  Suppose  that  f(x,u)  of  l  is  an  analytic  Devalued  function. 

I  is  locally  linearizable  at  (0,0)  if 

(i)  the  elements  of  E  are  linearly  independent, 

(ii)  span  E..  *  span  (E^HE),  for  1  < i  <m,  and 

(iii)  there  exists  k(  <»)  such  that  y^espan  ((F^)k)  for  1  ±1  £m  and 
i  >1,  where  (F^)k  is  composed  of  the  first  k  columns  of  D^. 

Given  the  system  (4)  we  choose  the  Kronecker  indices  {K^ in  a 
similar  way  to  the  continuous  time  case  (Hunt  and  Su  1981).  First  we  form 
the  matrix 


/  af  \ 

^u/(0,0) 

/af\  /  af  % 
'dx'(0,0)'dUj  (0,0) 


(JL\ 

lau2;(o,o) 


/  af  \ 

•**  v3um;(0,0) 


vax;(0,0)l3u2;(0,0)  •••  '3x'(0,0)'3um'(0,0) 


'3x'(0,0)'3u-j  '(0,0) 


/3_f\^  1  ,Jf_)  /_af\^  ^  (  df  \ 

'3x'(0,0)'3u2  (0,0)  '3x'(0,0)'3um'(0,0) 


Set  a.  =  number  of  linearly  independent  vectors  in  the  first  (i+1)  rows, 
for  0£i  <N-1.  Take  Yq  =Oq  and  y^  sc^  - _-j  for  1  £i  £N-1 ,  and  define 
k,.  to  be  the  number  of  Yj  with  y^  >i. 
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Appendix 

Lemma  A . 1 :  The  equation  (2)  holds. 
Proof:  Clearly, 
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where  (el  „}  are  to  be  determined.  Then 
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By  (A.l ) ,  since  bJ>0  *  D^,  B^q  =  for  m>1. 

true  when  m=l.  Now  suppose  that  (2)  is  true  for  m<p. 
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Note  that  (2)  is 
Let  1  <t  £p.  Then, 
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Changing  the  dummy  variables  k^ .kg,. . . ,k^  of  the  second  term  Into 
k2.k3,. . . ,k^,  respectively,  the  second  term  becomes 
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Thus  (2)  is  true  for  m=p+l  and  1  l^£p.  By  (A. 3),  it  is  easy  to  see  that 
(2)  is  true  for  m*£=p+l.  Hence  (2)  Is  true  for  m=p+l.  By  induction, 
(2)  is  true  for  m  >J . 
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Let  h(x):  DT-*  IR  be  a  C  function. 


Lemma  A. 2:  If 
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Proof:  (A. 4)  Is  equivalent  to 
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Therefore,  d‘  e  span  (Q) ,  where 
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Now  consider  the  matrix  defined  by 
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where  1  £a^  and  1  <a2<N.  Clearly  Image  ($p  *  span  (Q).  Similar 
arguments  can  be  applied  for  " ,(Dx+1h)x-0*  Therefore* 
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